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ABSTRACT 


Theoretical arguments for a 30% increase in the solar luminosity over 
the past 4.7 billion years are reviewed. A scaling argument shows that this 
Increase can be predicted without detailed numerical calculations. The 
magnitude of the Increase is independent of nuclear reaction rates, as long 
as conversion of hydrogen to helium provides the basic energy source of the 
Sun. 


The effect of the solar luminosity Increase on the terrestrial climate 
is briefly considered. It appears unlikely that an enhanced greenhouse effect, 
due to reduced gases (NH^, CH^), can account for the long-term paleodimatic 
trends. 



Research supported in part by KASA grant NAG 5-13 and NSF grant AST 79-19688. 


INTRODUCTION 


Climatically significant changaa of the solar luminosity (L) have been 
postulated to occur on time scales ranging from a few years to billions of 
years. The shorter time scales have been discussed extensively a: this 
conference. In the present review, I will restrict myself to the longest time 
scales (SJlO^ yr.) and discus* the basis for the astrophysical conclusion that 
the Sun was ~ 30% fainter 4.7 x 10 9 yr. ago and that the evolution since the 
Sun's formation requires a slow, but steady, increase in L. 

I should note that this is the only change in L predicted by stellar 
evolution theory, in its standard form. This prediction is common to all 
modern calculations and is supported by a large body of data from observational 
stellar astronomy (see reference 1 for a review of the observational evidence) . 
Nevertheless, the validity of this result has been questioned because of the 
apparent conflict with proxy Indicators of the Earth's past climate (ref. 2-4). 
For this reason, a review of the theoretical arguments for the long-term 
increase of L is in order. 

Stellar evolution is governed by nonlinear differential equations derived 
from conservation laws and considerations of energy transport processes. 
Analytic solutions do net exist for any cases relevant to the Sun so numerical 
solutions must be used. Modern calculations require complex computer codes 
Incorporating a variety of physical data on nuclear parameters, transport 
coefficients, and thermodynamic properties. In this respect, the situation is 
similar to that encountered in current theoretical investigations of the 
terrestrial climate. This may, in fact, explain the reluctance of the climato- 
logists to accept the astrophysical result: climatologists understand the 
pitialls of accepting solution < obtained from complex computer codes at face 
value. For this reason, I will largely avoid discussion of numerical models 
and base the astrophysical case on simple scaling laws. 


A SCALING MODEL OF THE SUN 

We begin by requiring that the Sun be in hydrostatic equilibrium, with 
gravitational forces balanced by the pressure gradient. The free-fall time 
of the Sun Is on the order of an hour and any departures from hydrostatic 
equilibrium would show up as luminosity and radius changes on this time scale. 
For the spherically symmetric case, hydrostatic equilibrium is expressed as 
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whara C is Che gravitational constant, P and p are the pressure and density at 
a distance r from the center and m is the mass interior to r. Measurements 
of the visible solar disk show thaE the Sun is spherical to Within 1 part in 
10 5 (ref. 5). 

We can construct a one-zone model by replacing (1) by a finite-difference 
equation evaluated between the center c and aurface a, with mean values 
enclosed in brackets ( ) : 
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Applying the boundary conditions P ■ 0, r ■ R (radius), and r ■ 0 gives 

8 8 C 

m p 


P - G 
c 


< + > *• 


( 2 ) 


The scaling laws for the mean values are: 
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where m is the total mass. Inserting these scaling laws Into equation (2) 
gives 
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To proceed further, we need an equation of state, relating the pressure 
to the density and temperature (T) . For typical conditions characterizing 
the bulk of the solar interior (p at 1 g/cm^ anc. T ea 10^ to 10^ °K) , Coulomb 
interaction energies are at least 2 orders-of-magnitude smaller than particle 
kinetic energies. Thus, the ideal gas law is an excellent approximation and 
this is what differentiates a star from a planet. Applying the ideal gas law 
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to the canter gives 
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where k in the Boltzmann constant and the particle density is expressed as the 
Avogadro number N. divided by the mean mass p (in atomic mass units) per free 
particle. Eliminating P between (6) and (7) and noting that the central 
density must scale as thS mean density gives 



If the scaling law (5) la used to replace R, we get 
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We now turn to the question of how energy is transported from the core, 
where nuclear reactions produce energy, to the surface. Due to the high 
temperatures, radiative transport of energy is very efficient and dominates 
over the bulk of the interior. The mean free path £ of a photon is typically 
1 cm so the photon diffusion approximation is valid to order £/R c* 10** . 

The radiative diffusion equation with spHevical symmetry is 

L - £T <„ 

r 3 kp dr 


where L is the total flux across a spherical surface at distance r from the 
center, r o is the btefan-Boltzmann constant and k Is the Rosseland -mean opacity 
coefficient. Again, we ;:?e the one-zone difference approximation to write 
this equation as 
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Application of our previous boundary conditions plus T^ ■ 0 (i.e. T^<< T^) 
givss 
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where, sines we are now dealing with a scaling law (proportionality' , 
can be replaced by L, r by R, etc. 

To evaluate the opacity coefficient *, we note that, in the solar 
interior, hydrogen and helium will be completely ionized and the heavier ions 
will be stripped of most of their electrons. Hydrogen and helium affect k 
through free-free transitions while the heavier elements contribute primarily 
through bound-free transitions. Both processes are reasonably repreaented by 
the hydrogenlc approximation so the absorption coefficient for a glvan ion 
varies inversely with the cube of the frequency. Although individual 
ionization states may contribute "noise" to the detailed dependence of k on p 
and T, the broad dependence i» given by Kramers' opacity: 
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Putting this result into equation (10) , and using equation (8) to eliminate 
the temperature, gives 


5.33 0.17 ,7.5 
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_^The present rate of mass loss, due to the solar wind, is roughly 
l<T XH m /yr. (ref. 6) and there is no reason to believe that the mass loss rate 
in the 0 past was great enough to significantly affect m. The density 
dependence in (12) is so weak that we may also neglect changes in this 
parameter. Thus, the luminosity is primarily dependent on the mean molecular 
weight y and we rewrite (12) as 



We let X, Y, and Z denote the fractional abundances, by mass, of hydrogen, 
helium, and metals (X+Y+Z ■ 1). For a fully ionized gas, 
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The mean molecular weight increases with time due to conversion of hydrogen 
(w “ 1/2) into helium (u ■ A/3), by nuclear reactions, producing Q ■ 6x10*® 
erg per gram of hydrogen consumed. This energy must supply the luminosity 
of the Sun. Since Xm is the total mass of hydrogen. 
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Differentiating (1A) with respect to time, and noting that 
dY/dt “ - dX/dt, we get 
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Finally, we can eliminate p between equations (13) and (16). The resulting 
differential equation can be directly integrated to give 
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Since nuclear reactions are confined to the core, the present photo- 
spheric abundances should reflect the initial composition. Thus, we may 
evaluate p using X =* 0.71 and Z - 0.02. Equation (17) becomes 
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where L is expressed in units of the present solar luminosity (taken as 
3.9x10^ erg/s) and t is the present solar age (A. 7x10^ yr.). The initial 
luminosity required tS match the present solar luminosity at t Q is L(o) ■ 0.76. 


6 



The scaling arguments predict that the Sun was initially 24% fainter than 
the present luminosity. A comparison of the luminosity evolution according 
to equation (18) with results from detailed numerical models (ref. 7) is 
shown in figure 1. As noted by D. 0. Gough, the evolution predicted by 
numerical models is accurately represented by 



e 


where L is the present solar luminosity. This formula, rather than equation 
(18), is recommended for studies of the evolution of planetary atmospheres. 


SUMMARY OF THE ASTROPHYSICAL CASE 


The above analysis shows that a quantitative prediction of the 
evolutionary increase of the Sut:'s luminosity may be made without detailed 
knowledge of the physical processes taking place in the Interior. Therefore, 
this prediction is not affected by the uncertainties in this knowledge. In 
particular, we did not have to specify any nuclear reaction rates since the 
net reaction rate, integrated over the solar mass, is determined by the 
measured solar luminosity. This is quite different from the case of the solar 
neutrino prediction, which is very sensitive to detailed nuclear reaction 
rates (ref. 8) . The discrepancy between the predicted and observed neutrino 
flux should not be used to argue that the luminosity prediction is also 
questionable . 


APPLICATION TO THE EARTH'S CLIMATE 


Sagan and Mullen (ref. 9) pointed out that an enhanced greenhouse effect 
due to higher concentrations of NH. and CH^ in the Earth's atmosphere, could 
have maintained a warm climate even with a tower solar luminosity. A similar 
conclusion was reached by Hart (ref. 10). This mechanism cannot, however, 
compensate for all of the solar luminosity evolution. 

F.ileological evidence (ref. 11) shows that the Earth's atmospheric 
chemist :y changed from reducing to oxidizing some 1.5 to 2 billion years ago 
and this would have removed the enhanced greenhouse effect due to reduced 
compounds. Roughly one*~half of the solai luminosity increase occurs during 
the last 2 billion years 1 " there is no evidence for a parallel increase in 
the Earth's mean surface w rature. Indeed, isotopic studies of Precambrian 
samples by Knauth and Epstein (ref. 12) indicate that the mean surface 
temperature has been d ecreasing during this time. Clearly, there is a need 
for further studies of the effects of crustal movements and volcanism. 
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biological activity, ate. on the long-term evolution of the Eerth'e climate. 
At present, it appears that the effects of solar evolution ia still buried 
in the "noise" due to other uncertainties in palaoclioatic models. 
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Figure 1 - Long-tern evolution of the tolar evolution. 

The evolution predicted by the eceling model {equation (18)] 
is ahown by the deehed line end the prediction from e detailed 
computer model ie shown by the eolid line. 





